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Abstract. The coherent evolution of the atom-cavity state in bimodal (cavity) 
experiments has been analyzed for a realistic time-dependence in detuning the 
atomic transition frequency. Apart from a 'smooth switch' of the atomic resonance 
from one to the second mode of a bimodal cavity, we considered also an 
additional (effective) interaction between the field modes of the cavity, known 
as 'communication channel'. Comparison of our model computations has been 
made especially with the measurements by Rauschenbeutel et al [2001 Phys. 
Rev. A 64 050301] who demonstrated for the first time the entanglement of 
the field modes in a bimodal cavity. It is shown that the agreement between 
the (theoretically) predicted and experimental phase shifts can be improved by 
allowing a 'communication' between the two field modes during a short but finite 
switch of the atomic transition frequency from one mode to the other. We 
therefore suggest that the details of the atom-cavity detuning should be taken 
into account for the future interpretation of bimodal cavity experiments. 
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1. Introduction 

During the last decades, entanglement has been recognized as a key feature of 
quantum mechanics that describes not only the correlation in composite quantum 
systems but is useful also for applications. Since the famous Bohr-Einstein debate 
and the seminal work of Einstein, Podolsky and Rosen [1] in 1935, indeed, a large 
number of entanglement studies helped improve our present understanding of nonlocal 
and nonclassical phenomena as they occur in the microscopic world. In quantum 
engineering and quantum information theory [5], moreover, entanglement has been 
found crucial for implementing new (quantum) information protocols, such as super- 
dense coding [3], quantum cryptography or even simple quantum algorithms [S]. 
However, despite of all progress in the design and description of entangled quantum 
system, their manipulation and controlled interaction with the environment still 
remains a great challenge for experiment owing to the fragile nature of most quantum 
states. Among various other implementations, an excellent control over the light fields 
and atoms has been achieved especially by using neutral atoms that are coupled to a 
high-finesse optical cavity [H [7] . 

In a recent experiment by Rauschenbeutel and coworkers [8], for example, the two 
modes of a superconducting (bimodal) cavity were prepared in a maximally entangled 
state by using circular Rydberg atoms. In this experiment, an entanglement of the two 
field modes was achieved, by properly adjusting the detuning of the atomic (transition) 
frequency of the Rydberg atom while it passes through the bimodal cavity. The 
entangled state produced in the cavity is probed later by a second atom that, after 
being detected, reveals the coherent evolution of the superposition of the cavity states. 
The coherence of the (two-qubit) cavity state was demonstrated by varying the delay 
time after which the probe atom interacts with the cavity. Different delay times give 
then rise to a oscillation in the final-state probability for the second atom to be found 
in either the ground or excited state of the two-level configuration, as supposed for the 
Rydberg atom. For this final-state probability, an analytical expression was derived 
by assuming an idealized time evolution of the quantum state of both, the photon 
field and the atoms, throughout the periods of atom-cavity interactions. Despite the 
high quality of the cavities today, however, such an ideal time evolution neglects a 
number of relevant effects, including the relaxation of the cavity field, the influence of 
external and internal stray fields, or imperfections due to the cavity mirrors and cavity 
geometry. In particular the proper matching of the atomic transition frequency to (the 
frequency of) one or the other cavity mode (the so-called atom-cavity detuning) plays 
an essential role on the superposition of the cavity states, because the atom has to 
mediate the interactions between the (cavity) modes while passing through the cavity. 

The experiments by Rauschenbeutel et al [8] nicely demonstrated how the field 
modes of a bimodal cavity become entangled or disentangled in a well controlled 
way by manipulating the (de-) tuning of the atom-cavity interaction. In contrast to 
the single-mode cavities, the use of bimodal cavities has been found an important 
step towards the manipulation of complex quantum states and for performing various 
fundamental tests in quantum theory [9l[T0l[Tll[T2l[13l[14l[15l[T6l[l2]. For these 
bimodal cavities, its important to know how (sensitive) the produced cavity state 
depends on the details and the particular shape of the detuning process. Since the 
two nearly degenerate frequencies of the cavity modes are fixed by the geometry of 
the cavity, a detuning of the atom (transition) frequency with regard to one of the 
field modes automatically affects also the detuning with regard to the other mode. In 



On role of the atom-cavity detuning in bimodal cavity experiments 



3 



our discussions below, we shall often refer to this detuning of the atomic transition 
frequency (of the Rydberg atoms passing through the cavity) with regard to the field 
modes briefly as the atom- cavity detuning. 

In the present work, we examine how the coherent evolution of the cavity state 
of bimodal cavities, i.e. the superposition of the two cavity modes, is affected by 
a realistic time-dependence of the atom-cavity detuning and how the atoms, when 
passing through the cavity, may interact with both cavity modes simultaneously. To 
exhibit these effects, we make use of two models for the atom-cavity detuning in which 
the (idealized) step-wise detuning of the atomic transition frequency from one to the 
other cavity mode is replaced by (i) a smooth detuning that happens in a short but 
finite ('switching') time interval and (ii) a simultaneous interaction of the atom with 
both cavity modes leading to a wave mixing in the cavity. We shall refer to these 
models as the (separate) single-mode interaction and communication channel model. 
In the latter model, the wave-mixing effects mainly arise due to imperfections of the 
cavity mirrors that causes the cavity modes to interact (communicate) with each 
other, of. Section 2. To follow the time evolution of the atom-cavity interaction, we 
combine analytical solutions for the Jaynes-Cummings Hamiltonian with numerical 
simulations, if the atom-cavity interaction is not resonant. For both models, the 
single-mode interaction and communication channel, the final-state probability for 
the 'probe' atom to be found in either the ground or excited state is then compared 
with the experiments in [8]. Using the cavity parameters from these experiments, 
we show that more realistic assumptions about the shape (and model) of the atom- 
cavity interaction improves the agreement between the theoretical predictions and 
experiment. A detailed account of the atom-cavity detuning will therefore play an 
crucial role also for manipulating and analyzing (more) complex cavity experiments 
in the future. 

In Section 2 we start from the Jaynes-Cummings Hamiltonian to describe the 
interaction of a two-level atom with the two field modes of a bimodal cavity. For 
this, we first (re-) derive the expression for the final-state probability in Subsection 
2.1 that the probe atom is detected in the excited state by using the matrix formalism 
and by assuming an idealized step-wise change in the atom-cavity interaction. This 
derivation sets the framework for introducing our more realistic models concerning 
the shape and explicit form of the atom-cavity interaction in sections 2.2 (single- mode 
interaction model) and 2.3 (communication-channel). In the latter model, we modify 
the Jaynes-Cummings Hamiltonian as to allow the atom to interact with both cavity 
modes simultaneously. In Section 3 then, the predictions from these models for the 
final-state probability are compared with the data of [8]. In particular, here we display 
and discuss how the 'switching time' affects the outcome of the experiments. Finally, 
a few conclusions are given in Section 4. 

2. Entanglement of two field modes using a bimodal cavity 

The use of the resonant atom-cavity interaction regime is the simplest way to generate 
entangled states between atoms and/or the field modes of a cavity. For a sufficiently 
high enough quality factor of the cavity mirrors, namely, this regime implies a 'strong' 
atom-field coupling for which the dissipation of field energy in course of the atom- 
cavity evolution becomes negligible. Such a small dissipation of the photon field plays 
a crucial role for engineering of coherent states in the framework of cavity QED. Apart 
from the quality of the mirrors, the correct matching of the atomic transition frequency 
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to the frequency of cavity field modes (the so-called detuning) is also an important 
ingredient in order to achieve the resonant atom-cavity interaction regime. 

In the following, let us first recall the basic components and notions of the cavity 
QED experiments by Haroche and coworkers [^. In these experiments, circular 
Rydberg atoms in quantum states with principal quantum numbers 50 and 51 are 
treated as 'two-level' systems, being in the ground state \g) or the exited state |e), 
respectively. Using the Stark shift technique, the frequency of atomic e g transition 
can be tuned in a well controlled way to the nearly degenerate frequencies of the two 
cavity field modes. If, for the moment, we consider only a single-mode cavity, then the 
evolution of the atom-cavity state is described (for both, a resonant and non-resonant 
interaction of the atom with the cavity) by the Jaynes-Cummings Hamiltonian [18j 



where loq is the atomic e g transition frequency, uji the frequency of the cavity field, 
and f2 the atom- field coupling frequency. In the Hamiltonian ([1}, moreover, oi and a'^ 
denote the annihilation and creation operators for a photon in the cavity, acting upon 
the Fock states \n), while S- and 5*+ are the spin lowering and raising operators that 
act upon the atomic states |e) and jg), which are the 'eigenstates' of the spin operator 
Sz with eigenvalues -1-1/2 and —1/2, respectively. Furthermore, if there is not more 
than one photon in the cavity, the overall atom-field state for a resonant atom-cavity 
interaction, i.e. for a zero detuning {Q ~ LOi^ ~ oji) , evolves according to [12] 



In the literature, this time evolution of the atom-cavity states is known also as Rabi 
rotation where t designates the effective atom-cavity interaction time in the laboratory 
and (n • t) the angle of rotation. Note that neither the state |e, 1) nor \g, 0) appears 
in the time evolution pa|) - (p6| . in line with our physical perception that the 'photon' 
of the e <-> (7 transition is 'stored' cither by the atom or the cavity, but cannot occur 
twice in the system. 

In contrast to single-mode cavities, a bimodal cavity possesses the feature of 
two nearly degenerate light modes of (usually) orthogonal polarization. Since the 
frequency of the cavity modes are fixed geometrically by the design of the cavity, its 
only the atomic frequency that can be tuned by means of external fields in order that 
the Rydberg atom interacts resonantly with either the first or the second cavity field 
mode. In practise, this detuning is done by applying a well adjusted time-varying 
electric field across the gap between the cavity mirrors [20] so that the desired (Stark) 
shift of the atomic transition frequency ujo{i) is achieved. By proper adjustment of 
ijjQit), one can induce the desired Rabi rotation of the atomic state interacting cither 
with regard to the first or the second cavity mode [5]. In fact, the development of 
bimodal cavities has been found important not only for the manipulation of complex 
quantum states but provides one also with an additional 'photonic qubit' in the 
framework of quantum information that may interact independently with the Rydberg 
atoms ('atomic qubits'). Below, we shall denote the two cavity modes by Mi and 
M2 and suppose that they are associated with the frequencies wi and W2, such that 
oji — UJ2 = 5 > Since, moreover, the frequencies of the two field modes are fixed, 
only the atomic e ^ g transition frequency can be changed and, as mentioned before, 




(1) 



|e,0) cos(17t/2)|e,0) -isin(m/2)|5,l), 
\g, 1) ^ cos {nt/2)\g, 1) - i sin {VLt/2)\e, 0) . 



(2 a) 
(26) 
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Figure 1. (a) Temporal sequence for producing and observing an entangled 
state of the two cavity modes. The pictograms are described in the text. The 
grey ellipse area indicates the place where the entangled state between the cavity 
modes is obtained. Plots (b) and (c) display the detailed shape of the atom-cavity 
detuning A(t) = uJo{t) - during the time periods [o, and [t, T-|- |g] , 
respectively. 



we shall briefly refer to the (de-)tuning of the atomic frequency with regard to the 
cavity modes as atom- cavity detuning. 

For a resonant interaction of the Rydberg atoms with both cavity modes, one 
need to be able to switch the atomic frequency from one to the other mode. Let 
A(t) denote the (time-dependent) atom-cavity detuning and let us assume that the 
atom is in resonance with the cavity mode Mi for A(i) = a;o(i) — uJi =0. It is 
then in resonance with M2 for A(t) = —5. For a sufficiently large frequency shift 
5 of the two cavity modes, moreover, the atom can be in resonance only with one 
of the modes and, hence, the time evolution of the (interacting) atom-cavity system 
A — Ml — M2 can be 'divided' into separate time evolutions owing to the A — Mi or 
A — M2 resonant interactions [H]. In practise, however, neither mode of the cavity can 
be frozen out completely because of the rather small shift S (« 2 in the experiments 
in Refs. [SlISlinDj the decoherence in the photon field [TUl as well as imperfections of 
the cavity mirrors, i.e. the local roughness and deviations from the spherical geometry. 
For any theoretical treatment of the decoherence and, especially, the imperfections of 
the cavity, we must proceed beyond the Jaynes-Cummings model ([T} by allowing the 
atom to couple with both cavity modes at the same time |1H 112] . Such a coupling 
then supports also the coupling between the cavity modes and is therefore referred 
to as communication channel in the literature [10], cf. Section 2.3. It is the goal of 
this work to demonstrate how a realistic atom-cavity interaction (coupling) can be 
analyzed quantitatively and can improve the theoretical predictions on the various 
probabilities, when compared with experiment [8] . 

2.1. The step-wise model for detuning the atom-cavity interaction 

With this short reminder on the development of bimodal cavities and the Jaynes- 
Cummings model for the atom-cavity interaction, we are now prepared to discuss all 
necessary steps for obtaining the coherent superposition of the cavity mode state as 
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reported recently in the experiments of Rauschenbeutel et al |3- In this reference, an 
idealized time evolution of the atom-cavity system was considered, and it was shown 
how the two cavity modes can be brought into the maximally entangled Bell state 

|vI/) = i=(eW'|0,l) + |l,0)) , (3) 

where |0), |1) denote the Fock states of the cavity mode Mi, and |0), |T) the Fock 
states of mode M2, respectively. 

In a step- wise model of the atom-cavity interaction, the coherent superposition ([3]) 
can be produced as follows. In the derivations below, we shall often apply the matrix 
representation of the time evolution operators in order to discuss the individual steps. 
The use of the matrix representation later facilitates also the application of our two 
models concerning a more realistic atom-cavity interaction. Suppose the cavity is 
initially 'empty', i.e. in the state |0, 0) = |0) x |0), and is crossed by the (so-called) 
source atom Ag prepared in the excited state |e). Being within the cavity, this atom 
is first tuned in resonance with the cavity mode Mi (A — 0) for a Rabi rotation 
Q,ti = 7r/2, and followed by a second rotation Q.t2 = tt, now being in resonance with 
the second mode M2 (A = —5). In the step- wise model of the atom-cavity interaction, 
it is essential that the detuning from mode Mi to M2 can be achieved instantaneously, 
i.e. within a period that is completely negligible to the time of interaction with the 
individual modes. This idealized sequence of atom-cavity interactions is shown in 
Figure [TJa) where the Rabi rotations are denoted by the black diamonds, containing 
the angle of rotation, and where the spatio-temporal evolution can be seen for both, 
the cavity modes as well as the atoms. Figure [IJb) displays the step-wise change in 
the atom-cavity detuning from the mode Mi to M2. After the time the source 
atom has passed the cavity (and turns out to be de-coupled from the state of the 
cavity modes). 

To understand the time evolution of the atom-cavity system as a whole, we can 
consider successively the resonant interaction of the atom with the two cavity modes 
Ml and M2. From Figure [T{b), we see that the (time-dependent) Hamiltonian can be 
written in the form 



, Hi{ujQ = LOi), < t < 7^ 

H.{t) = ei{t)Hl + e,{t)H2 -\rr ' ' (4) 

yH2{ujo-u2), 2n<^^2n 

with the two 'step functions' 



ei{t)^i-e(t-^, e^it) = 9{t 

and where {h = 1; /i = 1, 2) 

= LJoSz + ^ {S+af, + a+S-) + lo^ {a^a^ + ^ (5) 

refers to the Jaynes-Cummings Hamiltonian as discussed above. Following the text 
by Puri [19], here we introduce the operators 

such that [H^^N^ = 0. Especially the last relation enables us in the case that the 
atom is in resonance with the field mode [loq = oj^) to cast the Hamiltonian ^ into 
the equivalent but more convenient form 

= LJ^N^ + ^ (a+5_ + S+a^) . (7) 
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Of course, the Hilbert space associated to the Hamihonian (|3|) is given by the product 
space of the three (sub)systems: {\g), \e)), Mi (|0), |1)), and M2 (|0), For the 
sake of convenience, let us introduce an explicit notation for the 2^ basis states 



|Vi) = |e,0,0), IV2) = |g,l,T), 
IV5) = |e,l,0), IVe) = |g,0,l), 
which naturally fulfill the relations [i,j 



(V„|V, 



and 



IV3) = |e,l,l), 
IV7) = |e,0,l), 
,8) 



|v.,)(v, 



IV4; 



= 1.9,0,0), 
= |5,1,0), 



(8 a) 
(86) 

(9) 



Using these basis states, the wave function (of the atom-cavity system) can be always 
written in the form 



(10) 



and where Uij{t) = (Vi\U{t)\'Vj) refers to the matrix representation of the time 
evolution operator U{t) = cxp H{z) dz^. By this definition, the time-evolution 

matrix U^^\t) for the atom-cavity state a.t t = ^ as given by the time-dependent 
Hamiltonian ([4]) is factorized 



2n 



(- 



\2n) 



(11) 



where the operator U'^'^\t) refers to the evolution due to the As — Mi interaction and 
f7(2)(t) to that due to As - M2. Since, as mentioned above, the initial atom-cavity 
state is |Vi) = |e,0, 0), we find Ci(0) = Sn, and the factorization in pT|) is equivalent 
to the requirement 



/ ei{z)Hidz, f 02{z)H2 
Jo Ja 



dz 



= 0, 



(12) 



in accordance with the Hamiltonian (jl]). 

For the separate steps of the resonant atom-cavity interaction described by the 
Hamiltonian ([7]), the Schrodinger equation is given by solutions ([2a| - ((2&)) . Omitting 
the details of the derivation, for which we refer the reader to the literature [H] , the 
matrix representation of the operators U^^\t) and U'''^\t) in the basis (|8ap - ((86| is 
given by 
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(13) 
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and 



and where 



/ x'{t) 





v'{t) 



V 



x{t) = x{t) = COS ( ^ ) , 
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1 
1 / 



(14) 



y{t) = m 



-I sm 



(15) 

In this notation, the prime refers to the additional phase factor □'(t) = e"^* • □(t) 
that arises from the energy difference 1%5 of the two cavity modes and is accumulated 

and 



during the {A — M2) Rabi rotation [5]. Substituting the expressions (|lip . ([T 
(|14p into the wave function (|10p and by making use of the initial condition Ci(0) 
we then obtain at the time t 



1^ two non-zero coefficients 



C6 



I* 



2fi 



|1,0) |5) 



(16) 



(17) 



and, hence, the total wave function (up to an irrelevant phase factor) is 

This wave function is equivalent to those derived in [5] and is the same as displayed 
above ([3]). Apparently, the two cavity modes form now are maximally entangled state, 
while (the state of) the source atom is factorized out. 

Unfortunately, the cavity state is inaccessible for direct measurements. In order to 
'prove' the coherent superposition p?)) of the two cavity modes, another (probe) atom 
Ap has to be sent through the cavity after a time delay (2^ ~ t^*^ experiments 

in [H], this time delay was chosen between about .. 710 jis. The purpose of the probe 
atom is to 'read off' the state of the cavity modes and to copy this information upon 
its own state. Similar as the source atom, the probe atom Ap interacts with the cavity 
during the time interval [T, T + being in resonance with the cavity modes Mi 
and M2 as shown in Figure [U^c). Using the arguments from above, the whole time 
evolution of the atom-cavity system from zero up to the time {T + ^) is given by 
applying the evolution matrix 

Stt 



U,, T 



E 



u. 



(3) 
ml 



T 



Stt 
20 



U, 



Stt 

m 



(18) 



In this matrix, is given by pT|) and refers to the interaction of the cavity with 



the source atom, while the other matrices 11^:^^, C/^^^ and U^j^' describe the three 
subsequent steps: the free time evolution during the time interval 

diag{l,l,l,l,l,e^'*,l,l}. 



r(5) 



(19) 



which is diagonal in the basis 



and where a relative phase shift e**^* is 



accumulated due to the energy difference hS of the two cavity modes. Afterwards, 
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during the interval [T, T + . the further evolution of the atom-cavity system is 
driven by the Hamiltonian [of. Figure [TJc)] 



Hpit)^93{t)Hi+9iit)H2 = 



T <t<T 
T 



- < t < T + — 



(20) 



Since the atom-cavity interaction is as well described by the Jaynes-Cummings 
Hamiltonian 
and 



([7]), the matrices U^^ {t) and {t) coincides with the matrices (|T3| 

occurs. Including this phase 



141) . except that in u[f {t) the relative phase e*** 



factor, U-f\t) then becomes 
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(21) 



Having the complete time evolution (US]), we can easily determine the probability 
P(T) to find the probe atom in the exited state |e). after it has crossed the cavity and 
hit the detector D. This is the (final-state) probability that has been measured during 
the experiment [8] and that provides us with the information about the superposition 
of the cavity mode states ([T7|) . In practise, of course, the atom-cavity interaction is not 
ideal and has to be replaced by some realistic model for the atom-cavity interaction. 
For the idealized sequence of interaction, however, the (complete) time evolution owing 
to the matrix p8)) at the time t = ^ + T implies only the two non-zero coefficients 



Cl 



C6 



T 
■in 



T 



T 



-e 217 



1 + e^e 



iST 



iST 



2 / 2 

and, hence, the final-state probability to find Ap in the excited state is simply 

" _ 1 + cos (cj T + (j)) 



P{T) 



(22) 



Obviously, this probability oscillates between zero and one with the frequency lu = 6 
and with the phase (shift) cj) = tt (5/217. In the experiments by Rauschenbeutel et al 
[5] and Raimond et al [5], the shift between the cavity mode frequencies takes the 
value (5/27r = 128.3 kHz, while the atom-cavity coupling constant is VI/2tt = 47 kHz. 



2.2. Changing the atom- cavity detuning in a smooth way 

In the previous Subsection, it was shown how the (two) cavity modes become entangled 
with each other by detuning the atom-cavity interaction in a proper way. For a step- 
wise detuning A(t) as displayed in figures [Ijb) and (c), this leads to the maximally 
entangled state pT)) of the two cavity modes after the source atom has passed through 
the cavity, leaving the atom factorized out in a product state. This 'sudden' change in 
the detuning A(i) is however not physically feasible, and further analysis is required 
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(C) 



Klis) 



t(^s) 




Figure 2. Temporal sequence of the atom-cavity detuning if the interaction of 
the atom is changed smoothly from one to the other cavity mode. Plot (a) display 
the atom-cavity detuning A(f) = ajo(f) — oJi during the interval |^0, where 
the transition frequency of the source atom is detuned from the cavity mode 
Ml to the mode M2; this 'smooth change' is inflated in plot (b) for the interval 
[ '20^ ' ^'2n ] ' '^^^ plots (c) and (d) display the same but for the atom-cavity 
detuning A(t) = a;o(i) — ij^i during the time interval j^T, T-|- , i.e. when the 
cavity interacts with the probe atom. 



to understand how the evolution of cavity states is affected by a more reahstic time- 
dependence of the detuning. In this sections, therefore, we shall consider a smooth 
'switch' from the A^j^ — Mi to the A^/p — M2 resonant interaction within a finite time 
period, which we will denote by r. For instance, according to the experimental setup 
[8] this period is of the length r/fi < 1 fj,s, which corresponds to ^ 7 % of the overall 
1^ interaction time of the atom with the cavity, or about an angle jq in units of Rabi 
rotations, which is not anymore negligible (as we considered in the previous model). 
Below, we shall assume a smooth behavior for this finite switch as displayed in figures 
l^a) and (c), where the cavity is non-resonant to the source atom during the interval 
^] and with the probe atom during [T + T + '^^] , respectively. This 

smooth switch of the atomic resonance from one to the other cavity mode can be 
described mathematically by the Hamiltonian 

Hi, 0<t<^ 
HIit) = { „ , f^,~.?r , (23) 

where H\ and Hi are the resonant Jaynes-Cummings Hamiltonian ([7|) from above, 
and where (i) and (i) are given by 

(t) = S,A{t) +ujiNi + ^ (a+5_ + S+ai) , (24) 
H+{t) = S, [A{t) + S]+ UJ2N2 + ^ {4S- + S+a2) . (25) 
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Here, the superscript 'r' in (j23p is introduced in order to distinguish the Hamiltonian 
from that of the step- wise detuning (g]). Note that, for the present smooth change 
of the atom-cavity detuning, Hi and H_ [t) together describe the (time) evolution of 
the Ag—Mi interaction up to the time That is the resonant part of the atom-cavity 
interaction due to Hi is shortened by the time interval while the remaining time 
up to ^ follows the non-resonant Jaynes-Cummings Hamiltonian H-{t) in line with 
figure Ufa)- (b). Similarly, the As — M2 interaction is first non-resonant and modeled 
by the Hamiltonian H+{t) for the interval ^^^], and becomes later resonant with 
the cavity mode M2, as seen on figure [2Ic)-(d). Obviously, the Hamiltonian (^5)1 
can be understood also as an extension of the 'step-wise' model from Section 2.1 by 
introducing an additional short period for the non-resonant atom-cavity interaction. 
For the Hamiltonian (|23p. therefore, the time evolution matrix still factorizes and can 



be expressed at t = as 



m = E (l - Q (Q < (^«) 



where the matrices U^p (t) and U^J\t) arc the same as in (|T3|) and p^ . and where 

the matrix u[^\t) has still the form (|13p but with x{t), x{t), y{t), and y{t) now being 
solutions of the differential equations [15]: 

= ^-xit) + iA(i) y{t), = ^y{t) - \m x{t), (27) 

= + \A{t) xit), = ^x{t) - U{t) y(t) . (28) 

Analogue, the matrix ulp\t) has the form but with functions x{t)^ x{t), y{t), 
and y{t) that are obtained from (P7|) and (|28p by replacing A(t) A(t) + 5. 

No exact analytic solutions are known for the relations (P7|) - (P5| . Therefore, 
in order to derive the time evolution matrix (|26p . we need to solve these equations 
for x{t), x(t), y{t), and y{t) numerically up to < = and then to substitute their 

values into the matrices u'>f'\t)- As discussed above, the final-state probability to 
find the probe atom in the state |e) is obtained by performing a similar procedure for 
the Ap — Ml — M2 interaction during the time interval [T, T -\- |^] . During the time 
[III' ' course, the 'free' evolution of the cavity state still follows the trivial matrix 
(|19p . Putting all these pieces together, hence, the complete time evolution of the wave 
function for such a smooth change in the atom-cavity interaction is described by the 
matrix 

'3 V ^ 2Vi ) ^ 2n j^fe™V2i]J ™' \2n) \n 2nJ 

and where the matrices U-j\t), U-j\t), U^j\t), and ul'^\t) coincide with matrices 
(fT9|) . ((2T|) . mi), and ([261), respectively. For the initial condition Ci(0) = Sn, the 
final-state probability in model is again 

2 



p.(T) = |(Vi|v&,(r))|^ 



(30) 
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Figure 3. Temporal shape of the dimensionless functions fi(t) and /2(i) that 

/p — Ml interaction (solid line) and the 



describe the contribution of the A. 



M2 interaction (dot-dashed line) to the total Hamiltonian of the atom- 



cavity interaction. Plots (b) and (d) display also temporal shape of the product 
fl{t) • f2{t) (dashed line) that describes the contribution of 'communication 
channel' to the total Hamiltonian of the atom-cavity interaction. 



which, however, cannot be given in a closed form because of the need for solving the 
functions x(t), x(t), y{t), and y(t) numerically in the various matrices. Like the cosine 
function (j22p from the step-wise model, this probability still oscillates between zero 
and one with the frequency ujr and phase (j)^. In the Section 3, we display and discuss 
how the frequency and phase (shift) behave as function of the 'switching' parameter 
r, i.e. the time for which the cavity is non-resonant to the atom. 

2.3. Communication channel model: Allowing the mutual interaction of the cavity 
modes 

Up to the present, we have always assumed that, at a given time t, the atom interacts 
either with the cavity mode Mi or the mode M2, implying a sudden change in the 
interaction from one cavity mode to the other. In our notation above, this means 
that the overall atom-cavity interaction is understood as a set of pure Ag^p — Mi 
or Ag/p — M2 (non)resonant interactions. However notice that at the time ^ = ^ in 
figurelUb), for instance, the atomic transition frequency is equally far from the (nearly 
degenerate) frequencies of both cavity modes and, hence, we expect the atom to 'feel' 
the (non-resonant) contribution of Mi and A/2 simultaneously. As additional sources 
for such a simultaneous interaction also serve the imperfections of the mirrors, stray 
fields, decoherence, etc. Finally, this simultaneous interaction also implies a mutual 
interaction among the cavity modes to which one therefore refers as a communication 
channel 

To describe this 'communication' between the two cavity modes, we need to 
proceed beyond the Jaynes-Cummings model in order to allow the atom to couple 
with both cavity modes at the same time. In addition, in order to preserve the 
'correspondence' with the previous two models, this simultaneous interaction must 
happen only inside the time intervals [^^ff-i ^^^y-] and [T + ^2n^ ' ^ ~^ ^2n"^ ] ' gi'^^n 
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by the same 'switching' parameter t as in Section 2.2. Below, we shall model the 
total (source) atom-cavity interaction within the interval [O, by means of the 
time-dependent Hamiltonian 

H:it) = lHAt), ^<t<^ (31) 

I ^2, ^ < ^ ^ l§ 

with 

H^{t) = h{t)H^{t) + f2it)H+{t) + h{t)f2{t)Hi , (32) 

Hi = X {aia2 + ajoi) , (33) 

and where and H±{t) denote the Hamiltonians ([7]) and respectively. The 

temporal behavior of the detuning A(t) used in H±{t) is the same as we considered 
in the previous section and which is displayed in Figure [2l The two (dimensionless) 
functions /i(i) and f2(t) displayed in Figure a), thus determine the contribution 
of the non-resonant Jaynes-Cummings Hamiltonians H±(t) and the communication 
channel Hamiltonian Hj to the overall time-dependent Hamiltonian PT|) . In this 
model, again, the resonant As — Mi and Ag — M2 interactions are shortened by 
time interval r/fi and the evolution of the (atom-cavity) state is governed by the 
Hamiltonian (t) with the non-resonant part H± (t) and the communication channel 
Hj. The Hamiltonian Hj is the simplest expression that fulfills the above requirements 
and preserves the cavity field energy. The contribution of Hj to the Hamiltonian (t) 
is given by the time-dependent dimensionless product /i (t) /2 {t) . The shape of this 
coupling term is displayed in Figure [3Kb) by a dashed line, which has its maximal 
value 1/4 at tmax = 7r/2r2. Making use of the functions ffj,{t) in the Hamiltonian ([3T|) . 
we expect to describe the atom-cavity interaction during the switching period r/fi in 
a realistic fashion. Moreover, below we shall consider also the coupling strength A of 
the cavity modes such that the 'effective' coupling X fi{t) f2{t) satisfies the relation 

max max ) = n. (34) 

The meaning of the last relations can be understood as follows: The strongest M1 — M2 
effective coupling should be (or, at least, should not exceed) the composition of the 
As — Ml and As — M2 couplings given by the values of fl/2, respectively. A simple 
calculation implies the value of A to be of the order of ^ 1 MHz. For this value, indeed, 
a good agreement with the experimental data of [8] has been found within the expected 
time (r/n < 1 fis), i.e. for that time which is needed to 'switch' the atomic frequency 
between the two cavity modes (see Section 3). Therefore, the Hamiltonian PT|) exposes 
the right 'correspondence' property with regard to the models we have introduced in 
the previous two Sections. This can be readily seen if we assume fi{t) = Oi{t) and 
/2(i) = (^2{t)\ then, expression (pij) reduces to the Hamiltonian HJ{t) and, if 
we additionally impose the limit t ^ 0, the Hamiltonian (|3ip simplifies to Hs{t) in 
expression 

The time evolution matrix at ^ = |^ as described by the Hamiltonian ((3T|) is 
given by 

5r(i^)^|:^.?'(^5^)<^'(n)<(^)' ('^) 
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Figure 4. The relative frequency (a) and pliasc (b) of the final-state probability 
as obtained experimentally and displayed for the four time intervals 7^ in |8] . The 
values from the experiment are displayed by dots and are combined by a dashed 
line in order to guide the eyes of the reader. See text for further explanations. 



where U^j\t) and ul^\t) refer to the expressions and (fT4|) . and where the matrix 



elements of U^^ {t) must fulfiU the Schrodinger equation 



- J2 {y^\H>< {t)\Vk) Uifit) . (36) 



dt 

k 



This set of differential equations cannot be solved analytically and need to be 
integrated numerically during the interval t = ^. Indeed, a similar procedure has to 
be performed later also for the Ap — Mi — M2 interaction during the time [T, T + |^] 
as seen from figures |3|c),(d) and figures [2lc),(d). Combining the pieces together, the 
complete time evolution of the atom-cavity interaction therefore becomes 

(-)(t + —\- V ( (-\ c/(^) (- - —\ 



' 2nj \ 2n J \nJ \n 2nJ 

k,Tn,l,n 

^£'(^-1^) (i^ I. (37) 

where the matrices U-f\t), ul^\t), U^j\t), and U-j\t) coincide with the matrices 



P9p . (HJ), (O, and ([SS]) respectively. Moreover, the final-state probability is given 
for the initial conditions Ci(0) = Sn by 



PriT) = |(Vi|vI'.(T))| 



2 



(38) 



Although no closed expression can be given for this probability (due to the numerical 
parts in the time evolution) , it still oscillates between zero and one with the frequency 
ujr and the phase 0t- that is different from the two other models in Sections 2.1 and 2.2. 
In the next Section, we shall discuss how these frequency and phase (shifts) behave 
as functions of the 'switching' parameter r. 



3. Comparison with experiment and discussion 

The expression of the final-state probability describes the time evolution of an 
idealized experiment [8] with a step-wise change in the resonant interaction between 
the two cavity modes. This expression neglects of course a number of additional 
effects, such as the cavity field relaxation, the infiuence of external stray fields, or 
imperfections due to the cavity mirrors and cavity geometry. The most significant 
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distortion of the ideahzed probabihty in (|22|) is given by the cavity field relaxation 
during the interval i-e. during the 'free' evolution of the cavity field, since 

during the Ag — Mi — M2 and Ap — Mi — M2 interactions, the system evolves in the 
strong coupling (resonant) regime, and hence, the dissipation of the cavity energy is 
negligible. This field relaxation mainly arises from the interaction of the cavity with 
the environment and its effect on the final-state probability ((22|) has been investigated 
recently by de Magalhaes and Nemes [10]. In this reference, the environmental degrees 
of freedom were modeled by a infinite set of harmonic oscillators whose annihilation 
(creation) operators Cf, (c^) are coupled linearly to the Hamiltonian of the cavity field 
modes as follows 

< f < 321 




3tt 

T <t<T 



H{t) = { H^oupied, §i<t<T (39) 

20 



with 



00 



Hcoupied = ^ a'l^a^ + ^ (uJh c+ct + 71b [aic+ + a'^ct] + y2b [a2C^ + at^b] ) 



/ . "^M ""M ""M ' 

fj.=l,2 b 

and where ujb describe the frequencies and the respective coupling constants. 
Following [TU], the final-state probability can be cast into the form 

P(T) = ^(r)l+cos(^^T + 0) ^ ^^^^^ ^^^^ 

with 

^(r) = e-("+««, B{T)='- ^ , ^^T-g, (41) 

and in which a, (3 denote the dissipation constants of the cavity, while the frequency lo 
and phase (f> coincides with those of expression ((22|) . In fact, the probability (|40)) has 
been found in a good agreement with the experimental data of [8] with regard to the 
amplitude and frequency of the oscillations. However, less agreement was obtained 
for the phase shift of the final-state probability as given hy (j) = 716/20,. Note that, 
according to (|40p . the two arguments of the cosine: uj and (f>, are independent of the 
dissipation constants a and /? of the cavity. Therefore, we can describe the frequency 
and phase shifts in the final-state probability separately from the damping of the 
signal, as given by functions (|1T|) . For the two models ([50)1 and ([35]), indeed, we will 
show below that a finite 'switch' (Section 2.2) and the communication between the 
cavity modes (Section 2.3) mainly affects the phase of the amplitude ([^D]) . as have 
been prognosticated in [TUj. Therefore, combined with the 'damping functions' ([^T]) . 
our investigations enables us to further improve the agreement with experiment. 

We now compare the predictions for the frequency and phase of the final-state 
probability, obtained with the models from Section 2, with those from experiment 
[8]. In these measurements, the probability for detecting the probe atom Ap in the 
excited state |e) was recorded and displayed for the four time intervals Ii = [48; 57] fis, 
I2 = [200; 207] (is, I3 = [400,408] /is, and h = [699,706] fis, respectively. From the 
observed probability (cf. Figures 2(a)-(d) in [5]), the four 'experimental' frequencies 
and phases have been extracted by a fit to the analytical function ([^O]) . where the 
cosine arguments have been considered as unknown values, separately for each time 
interval {ly = 1, . . . ,4), and are displayed by the dots in figure [3J For the sake 
of convenience, moreover, these experimental data were normalized to the reference 
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Figure 5. Comparison of the relative phase ^ for a smooth change of the 
atom-cavity detuning (cf. Section 2.2) with experiment, (a) Resuhs are shown 
for the 'switching' parameter tQ~^ = 3.6 fis, for which the best agreement with 
the experimental data (dots) is found. The theoretical values (crosses) are joined 
by dashed linos to guide the eyes of the reader. Figure (b) displays the relative 
phase ^ as a function of t (curves C^). Every curve corresponds to a phase 
(pT as obtained for the time intervals The vertical dashed line indicates the r 
value from plot (a). 



frequency uj = 6 = 2Tr x 128.3 kHz and phase cj) = 7TS/2il ~ 4.29 of the 'ideahzed' 
probability As seen from figure 21 clear deviations of the extracted experimental 
phase 0e occur with regard to the reference phase cf). 

To understand how a smooth change in the atom-cavity detuning A(t) affects the 
frequency and phase of the predicted final-state probability, a number of steps have 
to be carried out for both models of Sections 2.2 and 2.3: (i) to evaluate numerically 
the final-state probability Pt{T) or Pt{T) for the time delay intervals lu from the 
experiment; (ii) to fit the numerical results from step (i) to the analytical function 
(|40p in order to extract the cosine arguments for every time interval; and (iii) to repeat 
the previous two steps for different values of the 'switching' parameter r. The (cosine) 
arguments obtained from this procedure can then be compared to the experimental 
dots from figure U in order to determine the parameter r that agrees better with the 
experiment. 

Figure [S] compares the relative phase (j)r/(t> for a smooth change of the atom-cavity 
detuning as modeled in Section 2.2. In Figure[n)Ja), results are shown for the 'switching' 
parameter r/fi = 3.6 yus, for which the best agreement with the experimental data 
(dots) is found. In addition, figure [5Ub) displays the relative phase ^ as a function of 
T for each time intervals I^. As seen from this Figure, our mixed analytic- numerical 
treatment of the atom-cavity state gives (as expected) the same final-state probability 
(|22p in the limit r ^ 0. Note that the differences between Wr and uj is about two 
orders of magnitude smaller than the differences in the phases (and are therefore 
not displayed here). This confirms that the frequency of the final-state probability 
is less sensitive to 'undesired' effects in the atom-cavity interaction than the phase 
shift. The best agreement between our predictions and experiment is found for a 
'switching time' r/fJ = r/f2 = 3.6 /is that is far outside of what is expected for the 
switch of the atomic transition frequency from one to the other cavity mode in the 
experiment [8]. Therefore, despite of its rather large effect on the phase-shift of the 
final-state probability, the finite switching time cannot explain the observations from 
the experiment as taken alone. 

In the communication channel model from Section 2.3, in contrast, the atom 
interacts for a short period with both cavity modes simultaneously, leading to an 
effective 'communication' between the cavity modes. Figure [UJa) displays again the 
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Figure 6. The same as figure [5] but for the communication-channel model 
in Section 2.3, (a) Results are shown again for the 'switching' parameter 

tS7~^ = 0.33 /js with the best agreement with the experimental data (dots). The 
theoretical values (crosses) are joined by dashed lines to guide the eyes of the 

reader. Figure (b) displays the relative phase ^ as a function of r (curves C^). 

Every curve corresponds to a phase <^r as obtained for the time intervals . The 
vertical dashed line indicates the t value from plot (a). 



relative phase (/>r/0 for that switching parameter, r/ll = 0.33 /.ts, for which the best 
agreement with the experimental data (dots) is found. This value nicely fits to the 
expected time that is needed to 'switch^ the atomic frequency between the two cavity 
modes. Moreover, the relative phase (/)r/0 as function of r is shown in figure (H^b) 
for each time intervals separately. We therefore conclude that a finite switch 
together with an effective communication of the two cavity modes, i.e. the model of 
the communication channel, is appropriate for describing the atom-cavity interaction 
when the atom passes through the cavity. 

Let us mention once more that numerical simulations were performed only 'while' 
the atomic frequency is switched between the cavity modes (ri7~^ < 1 \is\ both for 
the single-mode (Section 2.2) as well as the communication-channel model (Section 
2.3). For all other periods, the time evolution of the atom-cavity state was followed 
analytically (|15p owing to a purely resonant atom-cavity interaction. The slight 
deviation between the predicted and observed phases [Figures |6l[a) and EJa)] clearly 
indicate that there is no perfect agreement with the experimental data. However, no 
error bars has been given on the final-state probability in [8|, and hence this slight 
disagreement should be later on refined in accordance with error bars which occur in 
cavity QED experiments. 

4. Summary and outlook 

The coherent evolution of the atom-cavity state in bimodal (cavity) experiments 
has been analyzed for a realistic time-dependence of the detuning of the atomic 
transition frequency. Apart from a 'smooth switch' of the atomic resonance from 
one to the other cavity mode during a short but finite time interval (Section 2.2), 
we considered also an additional (effective) interaction, a so-called communication 
channel, between the two cavity modes (Section 2.3). While, in the former case, the 
atom interacts at any time with just one mode of the photon field, either resonantly 
or non-resonantly, a simultaneous interaction with both field modes is modeled by 
means of the communication channel. The outcome of the experiment is the final-state 
probability for that the probe atom is observed in the (excited) Rydberg state; this 
probability is analyzed by combining the analytical solution of the Jaynes-Cummings 
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Hamiltonian with numerical simulations for all those time intervals during which the 
atom-cavity interaction is tuned from one cavity mode to the other. 

Comparison of our model computations have been made especially with the 
recent measurements by Rauschcnbcutel and coworkers [5]. To this end, a series 
of computations have been carried out for the final-state probability as function 
of the delay time T and for different (realistic) choices of the 'switching time' r, 
the model parameter, during which the atom is non-resonant to the cavity modes. 
Making use of typical cavity parameters as reported in [8] , it is demonstrated that the 
agreement between the predicted and experimental phase can be improved by allowing 
a 'communication' between the two cavity modes, cf. Figures [5] and [6l Together with 
the recent study by de Magalhaes and Nemes [10], who investigated the damping of 
the probability amplitude due to the coupling of the cavity to the environment, the use 
of the communication channel from above brings the theoretical predictions in good 
agreement with experiment. For the interpretation of future bimodal experiments, 
therefore, it seems appropriate to take into account both, the decoherence as well as 
the details of the atom-cavity detuning dynamics. 

In fact, several proposals have been reported for performing bimodal cavity 
experiments [III [El [l3l [HI [H [HI [H] . In the proposal by Zubairy et al [13], for 
instance, a bimodal cavity is utilized to realize a quantum phase gate in which the 
two qubits are given by the cavity modes. Based on this gate, the authors suggested 
a scheme that enables one to implement Grover's search algorithm. Another fruitful 
branch of the bimodal cavity applications represent proposals [HI [TH [HI [H] , where 
the schemes for engineering of various entangled states between the atomic and/or 
photonic qubits have been reported. All these papers ignore, however, the effects of 
the 'non-instantaneous' detuning of the atom-cavity interaction from one cavity mode 
to the other, as well as the effective interaction among the cavity modes. Further 
investigation are therefore required to better understand how well these schemes may 
work in practise. Finally, we mention the papers [H] [H] where it has been proven 
that the coupling of both cavity modes to a common reservoir, as in line with Hcoupied 
((M)) . induces the tunneling of a field state from one cavity mode to another mode of 
the same cavity, and thus, opens a way to implement the environment assisted (short- 
distance) teleporting inside a bimodal cavity. We note that in order to follow the time 
evolution of such quantum systems embedded into a reservoir or under the external 
noise and to analyze different (entanglement or separability) measures, a 'quantum 
simulator' has been developed recently in our group [21j that can be utilized for such 
studies in the future. 
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